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Parity and time invariance violating electric dipole moment of 205 Tl is calculated using the rela- 
tivistic Hartree-Fock and configuration interaction methods and the many-body perturbation theory. 
Contributions from the interaction of the electron electric dipole moments with internal electric 
field and scalar-pseudoscalar electron-nucleon (T,P)-odd interaction are considered. The results 



are d(^ ub Tl) = -582(20)d e or d( 2Uj Tl) = -7.0(2) x 10~ is C 
surements are discussed. The results of similar calculations for 
d( 133 Cs) = 0.76(2) x W- 18 C SP 



e cm. Interpretation of the mea- 
133 Cs are d( 133 Cs) = 124(4)d e or 
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I. INTRODUCTION 



Recent very sensitive experiment performed in Seat- 
tle [l| puts very strong constrain on the electric dipole 
moment (EDM) of mercury. It now reads d( 199 Hg) = 
(0.49±1.29 st at±0.76 syst ) x 1(T 29 e cm, which is sevenfold 
improvement of the previous result of the same group. 
This renews the interest on the sources of atomic EDMs. 
In our previous paper 0] we calculated the EDM of mer- 
cury and other paramagnetic atoms due to nuclear Schiff 
moment, (T,P)-odd electron-nucleon interaction and in- 
teraction of the electron electric dipole moment (d e ) with 
nuclear magnetic field. The EDM of mercury due to nu- 
clear Schiff moment was also considered in a recent paper 
by Latha et al [3|. Other contributions include, e.g. inter- 
action of the electron electric dipole moments with inter- 
nal electric field and scalar-pseudoscalar electron-nucleon 
(T,P)-odd interaction. The latter two sources of atomic 
EDM are strongly suppressed in mercury due to zero to- 
tal electron momentum, J = 0. They give rise to EDMs 
of atoms with closed electron shells only in third order 
of the perturbation theory, when magnetic dipole hyper- 
fine interaction is also taken into account. The strongest 
constrain on the strength of these (T,P)-odd interactions 
came so far from the thallium experiment [3| (see also 
review Q for a detailed discussion). However, signif- 
icant advance in the accuracy of the measurements in 
mercury [l[ has changed the situation. Now the con- 
strain on the scalar-pseudoscalar electron-nucleon inter- 
action which comes from mercury EDM measurements 
is five times stronger than those from thallium measure- 
ments while the constrains on the electron EDM differ 
two times only All these results rely on atomic cal- 
culations which provide the link between atomic EDMs 
and the fundamental constants of the (T,P)-odd interac- 
tions. Due to significant progress in measurements it is 
important to revisit the calculations as well for the sake 
of improving their accuracy and reliability. 

The third-order calculations for mercury will be the 
subject of future work. In present paper we perform 
calculations of the thallium EDM caused by the elec- 
tron EDM and the scalar-pseudoscalar electron-nucleon 



(T,P)-odd interaction. The latter effect was considered 
by Martensson-Pendrill and Lindroth [l3[ and Sahoo et 
al The results differ almost two times, which is 

probably significantly larger than assumed uncertainty 
of both calculations. 

The EDM of thallium caused by electron EDM was 
considered by many authors % i, [3, E3, EH- The 
results show strong dependence on electron correlations 
and change significantly depending on how many correla- 
tion terms are included. The most complete calculations 
were performed by Liu and Kelly [TlT ] using the coupled 
cluster approach. The result is in relatively good agree- 
ment with the semiempirical estimations of Ref. [3] • 

All previous calculations of the thallium EDM treated 
the thallium atom as a system with one external elec- 
tron above closed shells. In present paper we consider it 
as a three valence electron system by including 6s elec- 
trons into valence space. We use the configuration in- 
teraction technique combined with the many-body per- 
turbation theory (the CI+MBPT [3, El method). We 
demonstrate that all instabilities of the results are due to 
strong correlations between external 6s and 6p electrons 
and using the configuration interaction technique to treat 
these correlations accurately leads to very stable results. 
We use exactly the same procedure for both (T,P)-odd 
operators which is another test of the consistency of the 
calculations. Our final result for the electron EDM is in 
excellent agreement with the most complete previous ab 
initio calculations by Liu and Kelly [Tl| , while our result 
for the scalar-pseudoscalar electron-nucleon (T,P)-odd 
interaction is closer to the result of Martensson-Pendrill 
and Lindroth [l| and differs significantly from Ref. [12] • 



II. METHOD OF CALCULATION 

The Hamiltonian of the scalar-pseudoscalar electron- 
nucleon (T,P)-odd interaction can be written as 



H SP = l ^=AC SP loToPN (r), 



(1) 



2 



where G is the Fermi constant, A = Z + N is nuclear 
mass number, C SP = [ZC^ P + NC% P )/A, and 7 are 
Dirac matrices. 

The Hamiltonian for the electron EDM interacting 
with internal atomic electric field E; nt can be written 
as E3] 



(2) 



Summation is over atomic electrons. 

Atomic EDM caused by any of the interactions (|1I2[) 
is given by 



\T>\M)(M\H IF \0) 
Eq — Em 



(3) 



where |0) is atomic ground state D = — eJ2i r i is the 
electric dipole operator and H TP is the (T,P)-odd oper- 
ator. 



A. The CI+MBPT method 

To calculate the EDM of thallium we consider it as a 
system with three valence electrons above closed shells 
and use the CI+MBPT method 0, Gj| for the valence 
electrons. The EDM of the atom in the CI+MBPT is 
given by the formula very similar to ([3]) by with slightly 
different meaning of the notations. First, the many elec- 
tron states |0), \M) are now three-electron states in the 
valence space. Second, the summation in the electric 
dipole operator D goes over valence electrons only while 
contribution from atomic core is taken into account by 
modifying the single-electron operator d: d — > d + SV COIC , 
where <5Vc 0ro is the correction to the electron core poten- 
tial caused by external field. Closed shell core does not 
contribute to the EDM in the second order due to zero 
total angular momentum. 

To perform the calculations we need to go through the 
following steps: (a) generate a complete set of single- 
electron states; (b) build an effective Hamiltonian in the 
valence space; (c) calculate core polarization; (d) perform 
summation as in ^ over a complete set of three-electron 
states. Let's consider these tasks in turn. 

We use the V N ~ 3 approximation as in Ref. [HI]. The 
calculations start from the relativistic Hartree-Fock pro- 
cedure for the triple ionized thallium ion. This gives us 
the states and potential V COIC = V N ~ 3 of the thallium 
core. We use the B-spline technique [19] to generate a 
complete set of single-electron states. These states are 
eigenstates of the Dirac operator with the electron po- 
tential V N ~ 3 . We use 50 B-splines of order 9 in a cavity 
of radius 40as . 

The effective CI+MBPT Hamiltonian for three valence 
electrons has the form 
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(4) 



where hi is the single-electron part of the relativistic 
Hamiltonian 

Ze 2 

hi = cap + - l)m e c 2 + V N ~ 3 + S X , (5) 

r 

and hi is the two-electron part of the Hamiltonian 



M r i> r 2) 



i" - r 2 



+ £2(7-1, r 2 ). 



(6) 



In these equations, a. and (3 are the Dirac matrices, V N 3 
is the Dirac-Hartree-Fock (DHF) potential of the closed- 
shell atomic core (N — 3 = 78, Z = 81), and S is the 
correlation operator. It represents terms in the Hamilto- 
nian arising due to virtual excitations from atomic core 
(see Ref. [3, [H| for details) . S = corresponds to the 
standard CI method. £1 is a single-electron operator. It 
represents a correlation interaction of a particular valence 
electron with the atomic core. S2 is a two-electron oper- 
ator. It represents screening of the Coulomb interaction 
between the two valence electrons by the core electrons. 
We calculate Si for s-electrons using the all-order tech- 
nique developed in Ref. [2Cj. Si for p and d electrons 
as well as S2 are calculated in the second order of the 
many-body perturbation theory using the B-spline basis 
set described above. We use 40 lowest B-spline states up 
to Z max = 5 to calculate S. 

The same B-spline states are used to construct three- 
electron states for valence electrons. We use 16 lowest 
states above the core up to / max = 2 for this purpose. 
The basis for the ground state is generated by allow- 
ing all possible single and double excitations from two 
initial configurations 6s 2 6p and 6s6p6d The basis for 
even states is generated by allowing all possible single 
and double excitations from three initial configurations, 
6s 2 7s, 6s 2 6d and 6s6p 2 . Variation of the basis size in- 
dicate that it is saturated with respect to ri max but not 
completely saturated with respect to i max - However, the 
contributions of the states with Z. max > 2 are small and 
can be neglected at required level of accuracy. 

The three-electron valence states are found by solving 
the eigenvalue problem, 



H cB * v = E v * v 



(7) 



using the standard CI techniques. Calculated and exper- 
imental energies of a few lowest-energy states of Tl are 
presented in Table |U One can see that the inclusion of 
S (CI+MBPT) leads to significant improvement of the 
agreement between theory and experiment. 

To calculate transition amplitudes we need to take into 
account the effect of core polarization by external field. 
This is done by means of the time-dependent relativistic 
Hartree-Fock method (see, e.g. Ref. j21|) which is equiv- 
alent to the random-phase approximation, so we will use 
the term RPA for short. The RPA equations for an ex- 
ternal field operator F 



{hi-e c )5^ c = -{F + 5V^-^, 



(8) 
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TABLE I: Three- 
citation energies 



electron removal energy (RE, a.u.) and ex- 
fcm - ) of thallium. 



State 




Theory 


Experiment" 






CI 


CI+MBPT 




RE 




-1.9177 


-2.0677 


-2.0722 




2po 

r l/2 








0.0 




2po 
r 3/2 


6345 


8049 


7793 


6s 2 7s 


Sl/2 


23023 


26810 


26478 


6s 2 7p 


2po 
r l/2 


30635 


34496 


34160 




2po 
r 3/2 


31541 


35507 


35161 


6s 2 6d 


E>3/2 


32313 


36553 


36118 




2 D5/2 


32363 


36624 


36200 


6s 2 8s 


Sl/2 


34893 


39037 


38746 


6s 2 8p 


2 p o 
r l/2 


37572 


41714 


41368 




2po 

r 3/2 


37936 


42122 


41741 


6s 2 7d 


D3/2 


38048 


42359 


42011 




2 D 5 / 2 


38074 


42395 


42049 


6s 2 9s 


Sl/2 


39563 


43728 


43166 


6s 2 8d 


D3/2 


40796 


45931 


44673 




2 4 D 5/2 


41593 


45971 


44693 


Qs6p 2 


Pl/2 


37195 


43545 


45220 




T 3/2 


40797 


48339 


49800 




P5/2 


44665 


52779 


53050 



"Reference \22 



are solved self-consistently for all states in atomic core 
in the same V N ~ 3 potential as for the DHF states. The 
operator F is either the electric dipole operator or the 
operator of the (T,P)-odd interaction, or any other op- 
erator (e.g. hyperfme interaction). The correction to 
the core potential SVP~ 3 is used to calculate transition 
amplitudes 



El vw = (VJF 



(9) 



Here ^ v and ty w are three-electron states found by solv- 
ing the CI equations ((7]). 

Calculated and experimental values of the electric 
dipole transition amplitudes and magnetic dipole hyper- 
fine structure (hfs) constants A for low states of thallium 
which are relevant to the calculation of the EDM are pre- 
sented in Table ITT1 Calculation of the hyperfme structure 
is a good way to test the wave function on short distances 
which is important for the matrix elements of weak in- 
teraction. The data in the Table show that the accuracy 
of the calculation of the El-transition amplitudes and 
hyperfme constants of s and P1/2 states is within few 
percent. 

Finally, the last task we must be able to do to calculate 
the EDM is to perform the summation over complete 
set of three-electron states. We use the Dalgarno-Lewis 
method [13] for this purpose In this method, a correction 
8fy v to the three-electron wave function of the ground 
state v is introduced and the EDM is expressed as 

d atom = 2(Sy v \F + 5Vp- 3 \y v ). (10) 
Here F is either the electric dipole operator or the op- 



TABLE II: El transition amplitudes and hfs constants A of 
some low states of 205 T1. 



States 



Calc. 



Experiment 



7si 



El transition amplitudes (a.u 
/a 1.73 1.81(2) 

6d 3/2 2.23 2.30(9) 

Hyperfme structure constants A (MHz 
21067 21311 
11417 12297 



6P1/2 
6P1/2 

6P1/2 

7S1/2 



Ref. [23, 24] 
Ref. [23, 24] 



Ref. [25] 
Ref. [26] 



TABLE III: Static scalar polarizability Qo of the thallium 
ground state (a.u.). 



«val 



Total 



Other 



-0.67 



44.50 



48.81 



49.2 a ,50.4" 



"Reference 
"Reference 



erator of the (T,P)-odd interaction. The correction 6^ v 
is found by solving the system of linear inhomogeneous 
equations 



(H ett - E v )6* v = — (G + 5Vg- 3 )* v 



(11) 



Here G is another operator from the pair d, H TP . If both 
operators F and G are the same the electric dipole op- 
erator d, then the expression similar to (fit)]) gives static 
polarizability of the atom. Table IIIII presents the results 
of the calculation of the static scalar polarizability ao of 
the thallium ground state. Here a cole is the contribution 
of the thallium core to the polarizability, <5a core is the cor- 
rection to the core polarizability due to Pauli principle 
which forbids excitations from the core to the occupied 
6s and 6p states, a va i is the contribution of the valence 
electrons to the polarizability. The final result is in good 
agreement with other CI+MBPT [28| and coupled clus- 
ter 1291 calculations. 



III. RESULTS 

The results of the calculations of the EDM of thallium 
in different approximations are presented in Table |IV] to- 
gether with earlier calculations. As it was pointed out 
in Ref. @, [1, Q thallium EDM is very sensitive to the 
strong correlations between 6s and 6p electrons. This 
interaction is treated pretty accurately in the configu- 
ration interaction technique used in present work. In 
contrast, all previous calculations treated thallium as a 
system with one external electron above closed shells. 
Therefore, present results are significantly more stable 
than earlier ab initio calculations. 

The main source of uncertainty for present calculations 
comes from the core- valence correlations. Most of the 
core-valence correlations are included via second-order 
correlation operator E. However, there are small con- 
tributions like higher-order correlations, correction to E 
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TABLE IV: EDM of Tl due to electron EDM (d e ) and scalar- 
pseudoscalar electron-nucleon (T,P)-odd interaction. 





10- 18 C SP e 


cm Comments 






this work 


-614 


-7.33 


single-configuration, no £ 


-537 


-6.43 


single-configuration with £ 


-625 


-7.49 


single-configuration 






in the ground state, with £ 


-602 


-7.22 


full CI but no £ 


-581 


-6.88 


full CI+MBPT but no RPA 


-582 


-6.98 


full scale calculations 


-582(20) 


-7 0(21 


final 






nf Vipt Ipii l?if ioti 


-585 




Rpf* N ll pnimlpiH plimtpr 


-1041 




Ref. [8], DHF+lst order MBPT 


-502 




Ref. [8], Tietz°+lst order MBPT 


-607 




Ref. [8], Green a +lst order MBPT 


-562 




Ref. [8], Norcross a +lst order MBPT 


700 




Ref. [6J, parametric potential 


-500 




Ref. [7J, semiempirical estimate 


-301 




Ref. [9J, 2nd order MBPT 


-179 




Ref. [10], 2nd order MBPT 




-4.056 


Ref. [12J, coupled cluster 




-7(2) 


Ref. [13], RPA+rescaling 



of correlations from Ref. [10] 



"Parametric potentials 

due to external field (structure radiation), renormaliza- 
tion of the wave function, etc. Quantum electrodynamic 
and Breit corrections are also expected to be small [30j . 
As one can see from Table IIVI the effect of including £ 



into full-scale CI calculations on the EDM of Tl is about 
3%. We use this as an estimate of the accuracy of our 
calculations. 

Similar calculations for cesium give the following re- 
sults (in agreement with previous calculations, see re- 
view [5]): 

d(Cs) = 0.759 x KT 18 C* SP e cm, (12) 

or 

d(Cs) = 124 d e . (13) 

The estimated error of these results is about 3%: 

The result of measurement of the EDM of 205 Tl [i] 
reads 

d( 205 Tl) = -(4.0 ± 4.3) x 10 _25 e cm. (14) 
Using the numbers from Table IIVI we find 

d e = (6.9 ± 7.4) x 10~ 28 e cm, (15) 

and 

C SP = (5.7 ±6.2) x 10~ 8 . (16) 

These numbers are in good agreement with the analysis 
of Ref. [|. 
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